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Abstract 

We present exact calculations of partition function Z of the q-state Potts model with next-nearest- 
neighbor spin-spin couplings, both for the ferromagnetic and antiferromagnetic case, for arbitrary tempera- 
ture, on n- vertex strip graphs of width Ly = 2 of the square lattice with free, cyclic, and Mobius longitudinal 
boundary conditions. The free energy is calculated exactly for the infinite-length limit of these strip graphs 
and the thermodynamics is discussed. Considering the full generalization to arbitrary complex q and tem- 
perature, we determine the singular locus B in the corresponding space, arising as the accumulation set 
of partition function zeros as n — > oo. 
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1 Introduction 



The q-statc Potts model has served as a valuable model for the study of phase transitions and critical 
phenomena |^. On a lattice, or, more generally, on a (connected) graph G, at temperature T, this model 
is defined by the partition function 

Z(G,g,z;)= ^e-/'^ (1.1) 

with the (zero-field) Hamiltonian 

^--^^'5..^, (1-2) 

where (Ji — l,...,q are the spin variables on each vertex i ^ G] (3 = {kBT)~^; and (ij) denotes pairs of 
adjacent vertices. The graph G = G{V, E) is defined by its vertex set V and its edge (=bond) set E] we 
denote the number of vertices of G as n = n{G) ~ \V\ and the number of edges of G as e(G) = \E\. We use 
the notation 

K = PJ , a = w~^=e-^, v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., v > corresponding to oo > T > for the Potts ferromagnet, 
and (ii) < a < 1, i.e., —l<v<0, corresponding to < T < oo for the Potts antiferromagnet. One defines 
the (reduced) free energy per site / = —f3F, where F is the actual free energy, via 

f{{G},q,v) = hm \n[Z{G,q,v)'/"] . (1.4) 

n — >oo 

where we use the symbol {G} to denote lim„^oo G for a given family of graphs. 

Let G' = (y, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex set V and an edge 
set E' C E. Then Z{G, q, v) can be written as the sum 

Z{G^q,v) = 

G'CG 

n{G) e(G) 
r=k{G) s=0 

where fc(G') denotes the number of connected components of G' and Zrs > 0. Since we only consider 



connected graphs G, we have fc(G) = 1. The formula (L5) enables one to generalize q from Z+ to M+ 
(keeping v in its physical range). This generalization is sometimes denoted the random cluster model j^]; 
here we shall use the term "Potts model" to include both positive integral q as in the original formulation 



in eqs. (1.1) and (1.2), and the generalization to real (or complex) g, via eq. (|l.l^). The formula (1.5) shows 
that Z{G^ q, v) is a polynomial in q and v (equivalently, a) with maximum and minimum degrees indicated in 
eq. (J^^)- The Potts model partition function on a graph G is essentially equivalent to the Tutte polynomial 
and Whitney rank polynomial |Q, ||], IT^-Q for this graph, as was discussed in [|l^ and is briefly 
noted in the appendix of this paper. 



In this paper we shall present exact calculations of the Potts model partition function for strips of the 
square lattice with next-nearest-neighbor (NNN) spin-spin interactions. Specifically, we consider strips with 
width Ly — 2 vertices and arbitrarily great length with various boundary conditions. To avoid increasing 
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the number of parameters, we take the nearest-neighbor and next-nearest-neighbor couphng strengths to be 
equal. This study is a natural continuation of our previous analogous calculations for strips of the square 
and triangular lattices H, n6[, and the reader is referred to these papers for background and further 



references (see also 18 ). We envision the strip as being formed by starting with a ladder graph, i.e. a 
2 X Ljc strip of the square lattice, and then adding an edge joining the lower left and upper right vertices of 
each square, and an edge joining the upper left and lower right vertices of each square (these two added edges 
do not intersect each other). Following our previous papers we shall denote this lattice as sqd, where 

the subscript d refers to the added diagonal bonds for each square. The Potts model with NNN spin-spin 
couplings on the square lattice and strips thereof can equivalently be considered as the Potts model with 
nearest- neighbor couplings on the sqd lattice with these diagonal bonds present. For our strip calculations, 
we take the longitudinal (transverse) direction on the strip to be the horizontal, x (vertical, y) direction, 
respectively. We use free transverse boundary conditions and consider free, periodic (= cyclic), and Mobius 
longitudinal boundary conditions. As we showed before |20[] , for a given value of — m, the Ly — 2 
cyclic sqd strip graph is identical to the corresponding strip with Mobius boundary conditions; hence, we 
shall refer to them both as Lm = sqd{{Ly = 2)m, FBCy, {T)PBCx), where the L here stands for "ladder" 
(with diagonal bonds added). The open strip will be denoted Sm- Following our labelling conventions in 
p2[ , Lx — m + 1 edges for an open strip Sm and = m edges for the cyclic/Mobius strip Lm- One has 
n-{Sm) — 2(m + 2), n{Lm) — 2m, e(Srn) = 6 + 5m, and e{L„i) — 5m. Each vertex on the cyclic/Mobius strip 
Lm has degree (coordination number) A = 5; this is also true of the interior vertices on the open strip Srm 
while the corner vertices have A = 3. Hence, the cyclic/Mobius strips Lm are A-regular graphs with A = 5, 
where a A-regular graph is defined as one in which each vertex has the same degree, A. For the infinite sqd 
lattice, A = 8. Note that, regarded as graphs, the cyclic/Mobius strips of the sqd lattice considered here, 
like the full 2D sqd lattice, are nonplanar (except for L,„ with m — 1, 2). In contrast, the open strips Sm of 
the sqd lattice are planar graphs. 

One interesting special case is provided by the zero-temperature Potts antiferromagnet. In general, for 
sufficiently large q, on a given lattice or graph G, the Potts antiferromagnet exhibits nonzero ground state 
entropy (without frustration). This is equivalent to a ground state degeneracy per site (vertex), W > 1, 
since Sq ~ kslnW. The T = (i.e., v = —1) partition function of the above-mentioned g-state Potts 
antiferromagnet (PAF) on G satisfies 

Z{G,q,-l) = PiG,q) (1.6) 

where P{G, q) is the chromatic polynomial (in the variable q) expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color ^, |Tl|, |23[ p4| . 
The minimum number of colors necessary for this coloring is the chromatic number of G, denoted x{G)- We 
have 

W{{G},q) = lim P(G,g)^/" . (1.7) 

n — yoo 

There are several motivations for the present study. Clearly, new exact calculations of Potts model 
partition functions are of value in their own right. A specific motivation is that this study provides exact 
results that reveal the effects of next-nearest-neighbor spin-spin interactions on the properties of the Potts 
model. It is of considerable physical interest what these effects are, since models with strictly nearest- 
neighbor interactions are only an approximation (albeit often a good one) to nature. For ferromagnetic spin- 
spin interactions (J > 0), the addition of (ferromagnetic) NNN interactions clearly enhances the tendency, 
at a given temperature, toward ferromagnetic ordering. For antiferromagnetic (AF) spin-spin interactions 
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(J < 0), the effect of the addition of (antiferromagnetic) NNN can be investigated by starting with the 
simple case of zero temperature. In the case 9 = 2, i.e., the Ising antiferromagnet, if one considers the square 
lattice at T = 0, there is complete antiferromagnetic long range order. However, in contrast, on the sqd 
lattice, the Ising model is frustrated. Closely related to this, on the sq^ lattice, the chromatic number is 
xi-^Qd) = 4, rather than the value 2 for the square (or any bipartite) lattice. 

For a regular lattice, as one increases the lattice coordination number, the ground state entropy of the 
g-state Potts antiferromagnet (if nonzero for the given value of g), decreases. This can be understood as a 
consequence of the fact that as one increases the lattice coordination number, one is increasing the constraints 
on the coloring of a given vertex subject to the condition that other vertices of the lattice adjacent to this one 
(i.e. connected with a bond of the lattice) have different colors. The addition of NNN spin-spin couplings 
to the Hamiltonian for the Potts antiferromagnet on the square lattice has a similar effect of increasing the 
constraints on the values that any given spin can take on, and hence decreasing the ground state entropy. 
Our exact calculations for the strips of the square lattice with NNN couplings give a quantitative measure 
of this effect. In a different but related direction, owing to the correspondence with the Tutte polynomial, 
our calculations yield several quantities of relevance to mathematical graph theory. 

Using the formula (1.5) for Z(G, q, v), one can generalize q from Z_|_ not just to R+ but to C and a from 
its physical ferromagnetic and antiferromagnetic ranges 1 < a < oo and 0<a<ltoaGC. A subset of 
the zeros of Z in the two-complex dimensional space defined by the pair of variables (g, a) can form an 
accumulation set in the n oo limit, denoted B, which is the continuous locus of points where the free 
energy is nonanalytic. This locus is determined as the solution to a certain {Gj-dependent equation in the 
two complex variables q and a [ p^ p^ . For a given value of a, one can consider this locus in the q plane, 
and we denote it as Bq{{G}, a). In the special case a — (i.e., v ~ —1) where the partition function is equal 
to the chromatic polynomial, the zeros in q are the chromatic zeros, and Sq({G},a = 0) is their continuous 



accumulation set in the n — > oo limit |25 -|2^. In a series of papers starting with |29(| we have given exact 
calculations of the chromatic polynomials and nonanalytic loci Bq for various families of graphs (for references 
on this a = special case, see H, ^ |20|). In particular, in we gave an exact determination of 
Bq for the [L^ oo limit of the) Ly = 2 sqd strip. A motivation for the present study is that it shows 
how the locus Bq that we calculated for the zero-temperature Potts antiferromagnet generalizes to finite 
temperature, as well as to the case of the Potts ferromagnet. With the exact Potts partition function for 
arbitrary temperature, one can study Bq for a ^ and, for a given value of g, one can study the continuous 
accumulation set of the zeros of Z{G,q,v) in the a plane; this will be denoted Ba{{G},q). It will often be 
convenient to consider the equivalent locus in the u — 1/a plane, namely Bu{{G},q). We shall sometimes 
write Bq{{G}, a) simply as Bq when {G} and a are clear from the context, and similarly with Ba and Bu- 
One gains a unified understanding of the separate loci Bq{{G}) for fixed a and Ba{{G}) for fixed q by relating 
these as different slices of the locus B in the space defined by {q, a) as we shall do here. 

Following the notation in |^ and our other earlier works on Bq{{G}) for a = 0, we denote the max- 
imal region in the complex q plane to which one can analytically continue the function W{{G},q) from 
physical values where there is nonzero ground state entropy as Ri . The maximal value of q where Bq 
intersects the (positive) real axis was labelled qc{{G}). Thus, region i?i includes the positive real axis for 
q > qc{{G}). Correspondingly, in our works on complex-temperature properties of spin models, we had 
labelled the complex-temperature extension (CTE) of the physical paramagnetic phase as (CTE)PM, which 
will simply be denoted PM here, the extension being understood, and similarly with ferromagnetic (FM) and 
antiferromagnetic (AFM); other complex-temperature phases, having no overlap with any physical phase. 
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were denoted Oj (for "other"), with j indexing the particular phase |32 . Here we shall continue to use this 
notation for the respective slices of B in the q and a or u planes. Another motivation for our present study 
is that it yields a deeper insight into the singular locus Bq for the zero-tenipcrature Potts antiferromagnet 
that we found for the ^ oo limit of the Ly — 2 sqd strip | |20[ | by showing how this locus changes as one 
increases the temperature from zero to finite values. 

We note some values of chromatic numbers for the strip graphs considered here: 



XiLm) 



XiSrn) = 4 

4 for even m > 2 

5 for odd m > 5 



(1.8) 
(1.9) 



Two degenerate cases are as follows: for ^ m = 2, the cyclic strip reduces to the complete graph|^ K^, 
while for TO = 3 it reduces to Kq, with xi^^p) — P- 

We record some special values of Z{G, q, v) below. First, 



Z{G,q = 0,t;) = 



(1.10) 



This implies that Z{G, q, v) has an overall factor of q and, in general (and for all the graphs considered here), 
this is the only overall factor that it has. We also have 



Z{G,q^l,v)^ V 
G'CG 



e(G') 



For temperature T — cxd, i.e., v — 0, 



Z{G,q,v = 0)^q'' 

rX(G)-i 



n 



U{G,q) 



ZiG,q,v = -l) = PiG,q)^ 

where U{G,q) is a polynomial in q of degree n{G) — xi^)- Hence, 

Z{G,q,v ^ -I) = P{G,q) ^0 for G = 5™, L,„ and (? = 1,2,3 



(1.11) 



(1.12) 



(1.13) 



(1.14) 



The result ( 1.14 ) implies that for q = 2 and 3, the partition functions Z{G, q = 2,v) and Z{G, q = 3,v) each 
contain at least one power of the factor (u + 1) = a; for g = 1, one already knows the form of Z{G, q = l,v) 
from ( 1.11 ). For the graphs with odd to, xi^m) = 5, so that Z{G,4:,v) contains at least one power of 
{v + 1) as a factor. 

Another basic property, evident from eq. (1.5), is that (i) the zeros of Z(G, q, v) in q for real v and hence 
also the continuous accumulation set Bq are invariant under the complex conjugation q q*; (ii) the zeros 
of Z{G, q, v) in v or equivalently a for real q and hence also the continuous accumulation set Ba are invariant 
under the complex conjugation a ~f a*. 

Just as the importance of noncommutative limits was shown in (eq. (1.9) of) on chromatic polyno- 
mials, so also one encounters an analogous noncommutativity here for the general partition function ( |l.5| ) of 
the Potts model for nonintegral q: at certain special points qs (typically qs = 0, l...,x(G)) one has 



lim lim Z{G, q, v 

n — ^oo q — >qs 



1/ri 



^ lim lim Z{G,q,v) 

q — >qs n — ^oo 



1/ra 



(1.15) 



^The complete graph Kp is the graph with p vertices each of which is adjacent to all of the other vertices. 
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Because of this noncommutativity, the formal definition (1^) is, in general, insufRcient to define the free 



energy / at these special points qg] it is necessary to specify the order of the limits that one uses in eq. 



(1.15). We denote the two definitions using different orders of limits as fnq and fqn- 



/„,({G},9,«)= lim \im n-HnZ{G,q,v) (1.16) 

n—>-oo q^Qs 



fqn{{G},q,v) = lim lim In Z(G, g, . (1.17) 

q — >qs n — >oc 



In Ref. [p9[ and our subsequent works on chromatic polynomials and the above-mentioned zero-temperature 
antiferromagnetic limit, it was convenient to use the ordering W{{G}^qs) ~ lim^^g^ lim„^oo -P(G, g)^/" 
since this avoids certain discontinuities in W that would be present with the opposite order of limits. In 
the present work on the full temperature-dependent Potts model partition function, we shall consider both 
orders of limits and comment on the differences where appropriate. Of course in discussions of the usual 



g-state Potts model (with positive integer q), one automatically uses the definition in eq. (1.1) with (|l.2|) 



and no issue of orders of limits arises, as it does in the Potts model with real q. As a consequence of the 



noncommutativity ( 1.1?: ) , it follows that for the special set of points q = qs one must distinguish between (i) 
{Ba{{G}, qs))nq, the contiuuous accumulation set of the zeros of Z(G,q,v) obtained by first setting q ~ qs 
and then taking n — > cxd, and (ii) {Ba{{G}, qs))qn, the continuous accumulation set of the zeros of Z[G, q, v) 
obtained by first taking n ^ oo, and then taking q ^ qg. For these special points, 

{Ba{{G}, qs)U + (Sa({G}, qs))qn . (1.18) 



From eq. (1.10), it follows that for any G, 

exp(/„,) = for g = (1.19) 

and thus 

[Ba)nq^% for g = 0. (1.20) 

However, for many families of graphs, including the circuit graph G„, and cyclic and Mobius strips of the 
square or triangular lattice, if we take n — > oo first and then g ^ 0, we find that {Bu)qn is nontrivial. 



Similarly, from ( 1.11 ) we have, for any G, 

(Sa)„, = for q=l (1.21) 

since all of the zeros of Z occur at the single discrete point a — (and in the case of a graph G with no edges, 
Z = 1 with no zeros). However, as the simple case of the circuit graph shows []l5| , {Bu)q7i is, in general, 
nontrivial. 

As derived in ||l^ , a general form for the Potts model partition function for the strip graphs considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated subunits (e.g. the 
columns of squares of height Ly vertices that are repeated times to form an x Ly strip of a regular 
lattice with some specified boundary conditions), is 

Z(G,g,«) =5]cg,,(Ag,,(<z,^;))" (1.22) 
i=i 

where N\ depends on G. This is a generalization of the result of |^ for the chromatic polynomial P(G, q) 
to the full Potts model partition function or equivalently, Tutte polynomial. 
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The Potts ferromagnet has a zero-temperature phase transition in the Lx ^ oo Umit of the strip graphs 
considered here, and this has the consequence that for cyclic and Mobius boundary conditions, B passes 
through the T — point u = 0. It follows that B is noncompact in the a plane. Hence, it is usually more 
convenient to study the slice of B in the u — 1/a plane rather than the a plane. Since a — s- cx) as T — > and 
Z diverges like this limit, we shall use the reduced partition function Zr defined by 

Zr{G, q, v) ^ a-''^^^^Z{G, q, v) = u''^^'>Z{G, q, v) (1.23) 

which has the finite limit Z^. — > 1 as T — > 0. For a general strip graph {Gs)m of type Gg and length — m, 
we can write 

Zr{{Gs)ra,q.v) = W^'C^^)-) ^ Cg„, (Ag. )" = ^ Cg„j ( Ag.,,,^)'" (1.24) 

with 

AG.,,,n=^i^«^=^'"^/'"AG.,, . (1.25) 



For the Ly = 2 strips of the sqd lattice of interest here, the prefactor on the right-hand side of ( 1.25 ) is 



2 Case of Free Longitudinal Boundary Conditions 

In this section we present the Potts model partition function Z(>5'm,g, u) for the Ly ^ 2 strip S„i with 
arbitrary length Lx = m + 1 (i.e., containing m + 1 squares) and free transverse and longitudinal boundary 
conditions. One convenient way to express the results is in terms of a generating function: 

oo 

T{S,q,v,z)=J2ziSm,q,v)z'^. (2.1) 

m— 

We have calculated this generating function using the deletion-contraction theorem for the corresponding 
Tutte polynomial T{Sm, x, y) and then expressing the result in terms of the variables q and v. We find 

T{S, q, V, z) = — — (2.2) 

V{b,q,v,z) 



where 
with 

and 



U{S,q,v,z)=As,o + As,iz (2.3) 

^5^0 = q{lQv^ + Ihqv^ + 15i;^ -f Av^q + Qvq^ + + + q^) (2.4) 
As.i = -qv'^{v + q){v + l){4qv + iv^q + 4w^ -|- 2v^ + 2q^ + Avq^ + v'^q^ + llgw^) (2.5) 

V{S,q,v,z) = I- (\2v^ + hqv + q^ + bv^ + l{)v^ + v'^)z 

+ v'^{v + l){2v^ + Av''^ + Aqv^ + v'^q^ + llqv^ +Avq + Avq^ + 2q^)z^ . (2.6) 
(The generating function for the Tutte polynomial T{S„i, x, y) is given in the appendix.) Writing 

2 

V{S,q,v,z)^\{{l-\s^,z) (2.7) 
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we have 



where 



As,(i,2) - 7; 



Tsi2 ± V Rsu 



(2.8) 



Tsi2 = 12u^ + 5qv + q^ + bv" + lOw'^ + v'' (2.9) 

and 

i?Si2 = 45v^ + 116v'' + 196v^ + lUv'^ + 224v^ + lQAqv^ + 4Qqv^ - Av^q^ 

+ ilv^q^ ~6v^q~10v^q + q^~2v^q'^ -lOv^q^ + 10vq^ + v^" + 10v^ (2.10) 

In p2| we presented a formula to obtain the chromatic polynomial for a recursive family of graphs in the 
form of powers of Xj 's starting from the generating function, and the generalization of this to the full Potts 
model partition function was given in [p^ . Using this, we have 

7{Q „ „.\ _ (^5,0^5,1 + ^S.l) , (^5,0^5,2 + Ag^i) 

^[Om,q,V) - — r r j H — r A5 2 [^■^^) 

(■^3,1 — ■^8,2) IA5_2 — Asaj 

(which is symmetric under As,i -^5,2)- Although both the Asj's and the coefficient functions involve the 
square root \/ Rsi2 and are not polynomials in q and ti, the theorem on symmetric functions of the roots of 



an algebraic equation ||33| , |34| guarantees that Z(Sm, q, v) is a polynomial in q and v (as it must be by (1.5) 
since the coefficients of the powers of z in the equation ( p.6[ ) defining these Agj 's are polynomials in these 
variables q and v). 

As will be shown below, in the limit m ^ 00 of this strip, the singular locus Bu consists of arcs that do 
not separate the u plane into different regions, so that the PM phase and its complex-temperature extension 
occupy all of this plane, except for these arcs. For physical temperature and positive integer q, the (reduced) 
free energy of the Potts model in the limit n — > cx) is given by 

/=^lnA5,i- (2.12) 

This is analytic for all finite temperature, for both the ferromagnetic and antiferromagnetic sign of the spin- 
spin coupling J. The internal energy and specific heat can be calculated in a straightforward manner from 



the free energy (2.12); since the resultant expressions are somewhat cumbersome, we do not list them here. 
In the r = Potts antiferromagnet limit v = —1, A5.1 ~ {q — 2)(q — 3) and A5.2 = 0, so that eq. ( |2.2| ) 
reduces to the generating function for the chromatic polynomial for this open strip 

rra 1 ^ g(g-l)(g-2)(g-3) 

T{S,q,v = -l;z) ^ — — — . 2.13 

1 - (9 - 2)(9- 3)z 

Equivalently, the chromatic polynomial is 

P(5,„,g) =q(q-l)[(g-2)(g-3)](™+i) . (2.14) 
For the ferromagnetic case with general q, in the low-temperature limit ?; ^ 00, 

As,i =v^ + 5v^ + 0{v^) As,2 = + 0{v) as v ^ 00 (2.15) 
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so that |A5,i| is never equal to |As,2| in this hmit, and hence Bu does not pass through the origin of the u 
plane for the n — > oo limit of this open sqd strip: 

u = O^Bu{{S}). (2.16) 

In contrast, as will be shown below, Bu does pass through u = for this strip with cyclic or Mobius boundary 
conditions. 



2.1 Bg{{S}) for fixed a 

We start with the value a — corresponding to the Potts antiferromagnet at zero temperature. In this case, 
Z{Sm, q,v — —1) — P{Sm, q), where this chromatic polynomial was given in eq. ( p. 14 ), and the continuous 



locus Bq = 0. For a in the finite-temperature antiferromagnetic range < a < 0.0790075, Bq consists of two 
self-conjugate arcs which cross the real q axis between 2 and 3. For 0.0790075 < a < 0.704942, Bq consists 
of two arcs which are complex conjugates of each other and do not cross the real q axis. The endpoints 
of the arcs occur at the branch points where the function Rsi2 in the square root is zero. As a increases 
through the value 0.704942, the two arcs pinch the real q axis with q smaller than 1, and Bq consists of a 
single self-conjugate arc crossing the positive real q axis and a short line segment on the q axis. As a reaches 
the infinite-temperature value 1, Bq shrinks to a point at the origin. In the ferromagnetic range a > 1, the 
self-conjugate arc crosses the negative real axis. In Figs. |l|, || and ^ we show Bq and associated zeros of 
Z in the q plane for the antiferromagnetic values a = 0.5 and a — 0.9, and the ferromagnetic value a — 2. 
Since the line segment in Fig. |^ is too small to be seen clearly, we note that it extends outward from the 
point at which Bq crosses the real axis and covers the interval 0.378027 < q < 0.379970. The line segment is 
sufficiently long to be visible in Fig. ^. 

2.2 B^{{S}) for fixed q 

We show several plots of the locus Bu for various values of q in Figs. H] - |[ In this section, for values 
of q where noncommutativity occurs, we display Bqn- Given the algebraic structure of Xsj, j = 1,2, the 
degeneracy of magnitudes \Xs,i\ = \Xs,2\ and hence the locus Bu occurs where (i) T512 = 0, (ii) Rsi2 — 0, 
(iii) if q is real, where Rsi2 < so that the square root is pure imaginary, and (iv) elsewhere for complex 
u, where the degeneracy condition is satisfied. The locus Bu does not enclose regions, and As,i is dominant 
everywhere in the u plane and degenerate in magnitude with As,2 on Bu- 

For large values of q, we find that Bu consists of two pairs of complex conjugate arcs, and a line segment 
on the negative u axis. The ten endpoints are the branch point zeros of \/Rsi2, and the line segment is the 
solution to the condition (iii) above. As q decreases, the locus B changes to consist of two self-conjugate 
arcs, as shown in Figs. || and |[ For 2 < q < 3, one self-conjugate arc crosses the positive u axis at the 
point where the condition (i) is satisfied. At g = 2, there is only one pair of complex conjugate arcs and a 
self-conjugate arc, the latter of which crosses the negative u axis. For 1 < q < 2, Bu consists of five disjoint 
arcs, as illustrated for the value q = 3/2 in Fig. ^. At q — 1, {Bu)qn is an oval (the solution to the degeneracy 
equation 2\u^{u — 1)^| = 1) that crosses the real axis at approximately 1.418 and —0.584. 

For < q < 1, B consists of two self-conjugate arcs, one pair of complex conjugate arcs, and a short line 
segment on the real u axis as illustrated for the value q = 1/2 in Fig. ^ (the line segment is too short to 
be clearly visible in this figure; it covers the interval 1.15474 < q < 1.15598). For q — 0, the locus {Bu)qn 
consists of two pairs of complex conjugate arcs. 
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Re(q) 



Figure 1: Locus Bq for the n ^ oo limit of the Ly = 2 sq^ strip {5} with free longitudinal boundary conditions, for 
a = ti + 1 = 0.5. Zeros of Z{Sm,,q,v = —0.5) for m = 20 (i.e., n = 44 vertices, so that Z is & polynomial of degree 44 in q) are 
shown for comparison. 
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Figure 2: Locus Bq-. same as Fig. |^ for a = 0.9 (i.e., v = -0. 
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Figure 4: Locus Bu for the n ^ oo limit of the Ly = 2 sqd strip, with free longitudinal boundary conditions, {S} with q = 10. 
Zeros of Z{Sm,q = 10,1") in u for m = 20 are shown for comparison. 
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Fi gure 6; Locus 3u'- same as in Fig. M for q = 3. 
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Fiff lire Ti Locus 3u'- same as in Fig. for q = 2. 
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Figure 8: Locus Bu'- same as in Fig. for q = 1.5. 
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Figure 9: Locus Bu'- same as in Fig. ^ for g = 0.5. 

3 Cyclic and Mobius Strips of the sqd Lattice 
3.1 Results for Z 

By using either an iterative application of the deletion-contraction theorem for Tutte polynomials and con- 
verting the result to Z, or a transfer matrix method, one can calculate the partition function for the cyclic 
strip graphs of arbitrary length, Z{G,q,v), G = Lm- We have used both methods as checks on the calcula- 
tion. Our results have the general form ( 1.22| ) with N\ = 5 and are 



A 



L,(2,3) 



— — T2?, ± \/ R 



23 



where 



T23 ^v'^ + 5v^ + Wv^ + Uv + 2q 
R23 = lUv^ + 32qv + Aq^ + + lOw^ + 45w^ + IQ&v'^ + ll&v"^ - I2qv^ - iqv^ + 22Av^ 



and 



Al,4 = As,i , Xl,5 — As,; 



where Xs,j, j — 1,2, were given above in (2 
The coefficients are 



CL,1 



9(9-3) 



CL.,2 = CL,3 ==9-1 
ClA = Cl,5 = 1 



(3.1.1) 
(3.1.2) 

(3.1.3) 
(3.1.4) 

(3.1.5) 

(3.1.6) 
(3.1.7) 
(3.1.8) 
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These can be expressed in terms of Chebyshev polynomials K 



.4) 

where Un{x) is the Chebyshev polynomial of the second kind, defined by 



c^'^^U,,{^^) (3.1.9) 



3=0 



J 



C^n(^) = E(-l)-'r /)(2^)"''' (3.1.10) 



where the notation [^] in the upper limit on the summand means the integral part of ^. The first few of 
the c''''''s are c'-^-' = 1, c'-^-' — q — 1, and c*^^^ — ~ 3q + 1. Thus, for the present case we have 



l(c(2)-c(°)) (3.1.11) 
cl.2 = cl_3 = c^^-*, and cl,4 = cl,5 = c^^'^ . The partition function Z[Lm, q,v) satisfies the general relations 



CL,1 2> 



(|1.10D-(|1.12D and (yj) 



Our main interest here is in large m and the m oo limit. However, for completeness, we make the 
following remark. If m > 3, then Lm is a (proper) graph, but the m = 1 and m = 2 cases requires special 
consideration; in these cases, Lm degenerates and is not a proper graph^. L2 is the multigraph obtained 
from the complete graph by doubling all edges except two edges that do not have a common vertex. Li 
is the pseudograph obtained by connecting two vertices with three edges and adding a loop to each vertex. 
Our calculation of Z{Lrm q,v) and the corresponding Tutte polynomial T{Lm, x,y) applies not just for the 
uniform cases m > 3 but also for the special cases m = 1, 2 if for m = 2 one includes the multiple edges and 



for m — 1 the multiple edges and loops in the evaluation of (1.1), ( |1.2| ), and (1.5). Note that in the T = 
case for the antiferromagnet, the resulting partition function, or equivalently, the chromatic polynomial, is 
not sensitive to multiple edges, i.e. is the same for a graph in which two vertices are connected by one edge 
or multiple edges; however, the general partition function (Tutte polynomial) is sensitive to multiple edges. 
The chromatic polynomial is sensitive to loops and vanishes identically when a pseudograph has any loops. 

3.2 Special values and expansions of A's 

We discuss some special cases here. First, for the zero-temperature Potts antiferromagnet, i.e. the case 



a — (v = —1), the partition function Z{L„nq,v) reduces, in accordance with the general result (1.6), to 
the chromatic polynomial P{Lm,q) calculated in In this special case 

Al,i = 2 (3.2.1) 

AL,2 = 2(3-g) (3.2.2) 

Al,3 = (3.2.3) 

XL,i ^ [q - 2){q ~ 3) (3.2.4) 

Al,5 = (3.2.5) 



(where these presume certain choices of branch cuts for square roots, so that, e.g., •\/(3 — q)^ = 3 — g, etc.). 

■^A proper graph has no multiple edges or loops, where a loop is an edge that connects a vertex to it self A multigraph may 
contain multiple edges, but no loops, while a pseudograph may contain both multiple edges and loops pil, f 
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For the infinite-temperature value a = 1, we have Xl.j = for j = 1,2,3,5, while \la = (f' , so that 



Z{Ljn,q,a = 1) 



^2m 



q", in accord with the general result ( 1.12 ). 



At g = 0, we find (with appropriate choices of branch cuts) 



A 



2v' 



{v + 'i){v^ + 2t; + 4) + Vl44 + + lOw^ + 45t;4 + 196^2 + I16u3 + 224u 



1 r 

2 



(u + 3)(w2 + 2t; + 4) - Vl44 + -y^ + lOw^ + 45u4 + 196w2 + 116t;3 + 224u 
Since there are two degenerate dominant terms, namely A2 = A4, it follows that 

q = is on Bq{{L\) V a . 



(3.2.6) 
(3.2.7) 
(3.2.8) 

(3.2.9) 



This was also true of the circuit graph and cyclic and Mobius square and triangular strips with Ly = 2 for 
which the general Potts model partition function (Tutte polynomial) was calculated in ||l5|. For q = 0, the 
coefficients c^^ ~ 0, cl^ = = —1, and cl^4 = cl,5 = 1 so that the equal terms cancel each other pairwise, 
yielding Z{Lm,q — 0,v) — 0, in accordance with the general result (1.10). The noncommutativity (1.15) 
occurs here: exp(/„^) = 0, while |exp(/q„)l = 

At g = 1,2 we again encounter noncommutativity in the calculation of the free energy: for q — 1, eq. 
( 1.11 ) yields the result Z{Lm, 9 = 1, a) = a^"\ whence 



exp(/„5({L},g- l,a)) =a^/2 

while fqn depends on which phase one is in for a given value of a. 
To discuss the special case q = 2, we first observe that 

lim lim Ar 3 lim lim Ar r 

Specifically, 

Al,3(« = -1) = if 9^2 
so that the left-hand side of ( ^.2.11 ) is zero, while 

Al,3(9 = 2) = 2«2 if 



(3.2.10) 

(3.2.11) 
(3.2.12) 

(3.2.13) 



so that the right-hand side of ( ^.2.11 ) is 2. Similarly, Xl,5{v = —1) = if g 7^ 2, but \L.5{q = 2) is nonzero 
if w 7^ —1 as given below. Thus, for the special case q = 2, with the understanding that in cases where the 
noncommutativity (|3.2.1l[ ) holds, we take the limit q = 2 first for general w, with v ^ —I, we have 



,2 



v{v + l){v + 2){v^ + 2v + 2) 



A 



i,(4,5) 



{v + l){v + 2) 



2v'^ + 2v + 2± \/v^ + 4t;5 + Sv'^ + 4v^ + iv'^ + 8v + i 



Further, for q = 2, cl,i = —1, while clj = 1 for 2 < j < 5; hence, the (Al^)™ and (Al^3) 
each other and make no contribution to Z, which reduces to 

Z{L„,,q^2,v)^ J2 (Al,,)" 

J=2,4,5 



(3.2.14) 
(3.2.15) 
(3.2.16) 
terms cancel 

(3.2.17) 
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Hence also, fqn fnq at q = 2. 

In order to study the zero-temperature critical point in the ferromagnetic case and also the properties 
of the complex-temperature phase diagram, we calculate the XL,j,uS corresponding to the A^j's, using eq. 



(1.25). In the vicinity of the point m = we have 

XL,i,u^2u^{l~-uf (3.2.18) 

and the Taylor series expansions 

Ai,2,n = l-u^ + 2{q- 2)u^ + O(u^) (3.2.19) 

Al,3,« = + 2{q - 4)w^ + 0{u^) (3.2.20) 

Al,4,„ = l + {q- 1)^" + 4(g - l)u^ + O(w^) (3.2.21) 

Al,5,« = + 4((7 - 2)u'^ + 0{u^) . (3.2.22) 

Therefore, at m = 0, Xl,2.u and Xl,4:.u are dominant and the boundary Bu is determined as the solution to 

the degeneracy of magnitudes of these dominant Alj^„'s, i.e., |Al_2,u| = |Al,4,„|, so that the point u = is 



on Bu for any q, with the understanding that for the values of q where the noncommutativity (1.15) and 
(1.18) occurs, we are referring to Bqn rather than Bnq. 



To determine the angles at which the branches of Bu cross each other at w = 0, we write u in polar coordi- 
nates as M = re**, expand the degeneracy equation |Ai,2,ti| = |Al,4.u|, for small r, and obtain qr^ cos(40) — 0, 
which implies that (for q ^ 0) in the limit as r = \u\ — s- 0, 

^^(2j+2V^ J -0,1,..., 7 (3.2.23) 

or equivalently, 6 — ±7r/8, ±37r/8, ±57r/8, and ±77r/8. Hence there are eight curves forming four branches 
on Bu intersecting at u — 0, with an angle of 7r/4 between each adjacent pair of curves at u = 0. The point 
w = is thus a multiple point on the algebraic curve Bu, in the technical terminology of algebraic geometry 
(i.e., a point where several branches of an algebraic curve cross [^). In the vicinity of the origin, u = 0, these 
branches define eight corresponding complex-temperature phases: the paramagnetic (PM) phase for — 7r/8 < 

< tt/8, together with seven Oj phases extending outward in the wedges {2j — l)7r/8 < < (2j + l)7r/8 
for j = 1, ..,7. These form one self-conjugate phase, 04 = 04, and the complex conjugate pairs of phases 

01 =0?, 02 = 01 and O3 = 05*. 

For q = 2 and for q = A the Potts antiferromagnet on the infinite-length, width Ly = 2 strip of the sqj, 
lattice has a zero-temperature critical point. In the q = 2 Ising case, this involves frustration. In order 
to study the T = critical point for the Ly = 2 strip for these two values of q, it is useful to calculate 
expansions of the Alj 's. Only Al,4 and Xl.2 are necessary for physical thermodynamic properties, while the 
full set of Xlj, j = 1, 2, .., 5 is, in general, necessary for the study of the singular locus B. 

For q = 2, besides the exact expressions Al,i = A^.a = 2{a — 1)^ and Al,2 = + a^, we have the 
expansions 

Xla = a + Aa^ + 9a^ + 0{a^) (3.2.24) 

Al,5 = 2a2 - 40^^ - 8a^ + O(a^) (3.2.25) 

As shown above, for fnq and Bnq, where one sets q = 2 first and then takes n — > 00, Xlj, j = 1,3, make 
no contribution, and Bnq is determined by the degeneracy in magnitude of the dominant terms among Xlj, 



16 



j = 2,4,5. From the expansions ( 3.2.24 ) and ( 3.2.25 ), it follows that in the neighborhood of the point 
a = 0, {Ba)nq is determined by the equation \Xl,2\ = \^l,4\- Writing a = pe*"^ and expanding the degeneracy 
equation |Al,2| = \^l,4\ for small r, we obtain 8p^ cos(30) — 0, which imphes that in the hmit as r = \a\ — > 0, 

^^(2j+2h^ J = 0,1,..., 5 (3.2.26) 

or equivalently, ((> = ±7r/6, (/> = ±7r/2, and ~ ±57r/6. Hence there are six curves forming three branches 
of Ba intersecting at a = and the angles between successive branches as they cross at this point are 7r/3. 
The point a = is thus a multiple point on B^- In the vicinity of the origin, a = 0, these branches define six 
complex-temperature phases: the paramagnetic (PM) phase for —tt/6 < < tt/6, together with the phases 
Oj for 1 < j < 5, with Oj occupying the sector (2j - l)7r/6 < 6* < (2j + l)7r/6. Note that O3 = O*, O4 = O*, 
and O5 = Ol- Xla is dominant in the PM phase and in the O2 and O2 phases, while Xl,2 is dominant in 
the Oi, 01, and O3 phases. 

For q — 4:, besides the exact expression Al.i = 2(a — 1)^, we calculate the expansions 

Al.2 = -2-2a + 3a2 + Cl(a3) (3.2.27) 

Xl,3 = a - + 0{a^) (3.2.28) 

Ai,4 = 2 + 14a-31a2 + 0(a3) (3.2.29) 

Al.5 = -3a + 33a2 + O(a^) . (3.2.30) 

It follows that in the neighborhood of the point a = 0, {Ba)nq is determined by the equation |Al_i| — |Al^4| 
and passes through a = vertically on the imaginary axis. 

3.3 Bg{{L}) for fixed a 

3.3.1 Antiferromagnetic case < a < 1 

For a — 0, i.e., the T — limit of the Potts antiferromagnet, the locus Bq consists of the union of the two 
circles 

: |(7-2| = 2 U |q-3| = 1 for a = (3.3.1) 

so that there Bq crosses the real q axis at g = 0,2,4, so that qc = 4; furthermore, the point g = 4 is a 
multiple point on Bq, where the smaller and larger circles intersect (shown in Fig. 4 in |20 ). This multiple 



point is a tacnode, i.e. the curves that intersect at this point have the same (vertical) slopes. In region Ri 
forming the exterior of the larger circle, |q — 2| > 2, Ai_4 is dominant. In region R2 forming the interior of 
the smaller circle, [g — 3| < 1, Al^i is dominant, while in region R3 comprised of the crescent-shaped area 
inside the larger circle and outside the smaller circle, i.e., with |g — 2| < 2 and |g — 3| > 1, Xl,2 is dominant. 

As the temperature increases from to infinity for the antiferromagnet, i.e., as a increases from to 
1, Bq contracts in to the origin, q = 0. Using our exact calculation of the Potts partition function for 
arbitrary T, we have determined Bq for general a. In Figs. [ic|-p^ we show some illustrative plots for the 
Potts antiferromagnet. As a increases from 0, the tacnodal multiple point that existed on Bq for a = 
disappears and the locus has the appearance illustrated by Fig. This is qualitatively similar to what we 
found for Bq for the {L^ — > 00 limit of the) Ly — 2 strip of the triangular lattice jl^ as a increased from 0. 
The middle crossing point at q = 2 remains fixed, independent of a, for a in the range < a < a„j, where 

am ^ 0.4294445 (3.3.2) 
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is the real solution of the equation 

2a^ + 3a^ + 3a-2 = (3.3.3) 

The reason for this is that this point is determined by the degeneracy equation |Al^i| = \^l.2\ where these 
are leading terms, and they cease to be leading at g = 2 as the right boundary sweeps leftward past this 



point; in turn, this occurs as Qc — 2, which yields the equation (3.3.3). As a increases from to a,„, region 
i?2 contracts in size, and finally disappears altogether as a increases through the value a,„. For a > a„n Al,4 
is dominant in the neighborhood of g = 2 as well as for q > 2. Again, this feature, that the crossing at g = 2 
remains fixed as a increases until the right-most part of B sweeps past it, thereby removing the region i?2, 
is qualitatively the same as what we found for the Ly = 2 strip of the triangular lattice jl^ . The value of 
a at which i?2 disappears in that case was (—3 + \/T7)/4 = 0.280776.... Note that B for the Ly — 2 square 
lattice strip is different in that in the a — case, qc — 2 and there are only two regions that contain intervals 
on the real axis: i?i for q < and q > qc ~ 2, and i?2 for < q < 2; since there is no third region including 
a line segment along the real axis, this, of course, precludes the phenomenon of the disappearance of such a 
region as a increases |l^ . 

The right-most point at which Bq crosses the real axis is given by the solution of the equation of degeneracy 
of leading terms |Al,i| = \^L,i\ for < a < am, 

4(1 -a)(a2 + 2a + 2) , 
(a + l)(a + 2) 

and is given by the solution of the equation |Al,2| — \^la \ for « > Om- The value of gc decreases monotonically 
from 4 to as a increases from to 1. For comparison, for the ^ oo limits of the Ly = 2 cyclic or Mobius 
strips of the square and triangular lattices, we found 0, ^ ^ 

qc{sq, 2 X cx), eye) = (1 - a){2 + a) (3.3.5) 

for all a, and 

qc{tm,2 X oo, eye) — (3.3.6) 

1 + a 

for < a < (1/4) (—3 + -^17) — 0.280776 (with a more complicated form holding for larger a). These all 
have the same property of monotonically decreasing to as a increases from to 1. 

3.3.2 Ferromagnetic range a > 1 

For the Potts ferroniagnet, as T decreases from infinity, i.e. a increases above 1, the locus Bq forms a lima- 



bean shaped curve shown for a typical value, a = 2, in Fig. 14, Besides the generally present crossing at 
q = 0, the point qc{{L}) at which Bq crosses the real q axis now occurs at negative q values. As was true 
of the model on the width Ly = 2 cyclic and Mobius strips of the square and triangular lattice jTB], |l^ , for 
physical temperatures, the locus Bq for the Potts fcrromagnet does not cross the positive real q axis. Note 
that this locus does have some support in the Re{q) > half plane, away from the real axis, which was also 
true of the analogous loci for the Ly = 2 cyclic and Mobius square and triangular strips. Finally, one could 
discuss the complex-temperature range a < 0; however, for the sake of brevity, we shall not do this here. 

3.4 Bu{{L}) for Fixed q 

We next proceed to the slices of B in the plane defined by the temperature Boltzmann variable u, for given 
values of q, starting with large q. In the hmit g oo, the locus Bu{{L}) is reduced to 0. This follows 
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Fi gure lOl Locus Bq{{Ly) for the n — ^ oo limit of the Ly = 2 sq^^ strip for a = 0.1. 
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Figure 11: Locus Bqi{L}): same as Fig. ^ for a = 0.25. 
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Figure 12: Locus Bq({L}): same as Fig. |lo| for a = 0.4. 
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Figure 13: Locus Bq({L}): same as Fig. |lo| for a = 0.5. 
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Fi gure 14; Locus Sq for the n — ^ oo limit of the sq^ strip with a = 2. 



because for large g, there is only a single dominant term, namely 

Al,4 ~ + hqv + 0(1) as q cxd . 



(3.4.1) 



Note that in this case, one gets the same result whether one takes q ^ oo first and then n = 2m — > oo, or 
n — > oo and then q —^ oo, so that these limits commute as regards the determination of S„. 

We first consider values of q 7^ 0, 1, 2, so that no noncommutativity occurs, and {Bu)nq = {Bu)qn = ^u- 
As discussed above, it is convenient to use the u plane since Bu is compact in this plane, except for the 
cases q — 2, and q — 4:, whereas Bu is noncompact because of the antiferromagnetic zero-temperature critical 
point at a = 1/u = 0. Extending the discussion in ||l^ to the case of the strip of the sqd lattice, we observe 
that the property that the singular locus Bu passes through the T = point u = for the Potts model with 
PBCx but not with FBCx, i.e., with periodic, but not free, longitudinal boundary conditions, means that 
the use of PBC^ yields a singular locus that manifestly incorporates the zero-temperature critical point, 
while this is not manifest in Bu when calculated using FBCx- 

For q = 10, the locus Bu is shown in Fig. |l^. Eight curves forming four branches on Bu run into the 
origin, w = 0, at the angles given in general in eq. (3.2.23). The Alj's that are dominant in these phases are 
Al,2 and Al,4, in an alternating manner as one makes a circuit around the origin. The PM phase includes 
the positive real u axis and evidently extends all the way around the curves forming Bu- The locus Bu also 
includes a line segment on the negative real u axis along which Xl,4 and Xl,5 are dominant and are equal in 
magnitude as complex conjugates of each other. In the two closed, complex-temperature O phases at the ends 
of the line segment, the dominant term is Al,2- As is evident in Fig. although the complex-temperature 
(Fisher p7|) zeros computed for a long finite strip lie reasonably close to the asymptotic Lx ^ 00 curve Bu 
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in general, these zeros have lower density on the curves running into the origin. 

For the q = 2 Ising case, the locus {Bu)nq is shown in Fig. |6|. One sees that, in addition to the eight 
curves intersecting at the ferromagnetic zero-temperature critical point m = 0, there is evidently another O 
phase that includes the negative real axis for u < — 1, and two complex conjugate pairs of O phases extending 
toward the upper and lower left, and upper and lower right. The complex-temperature phases in the vicinity 
of u = were determined above after eq. (3.2.23). Since the Ising antiferromagnet on the ^ oo limit of 
this strip has a zero-temperature critical point, it is useful to display the singular locus Ba in the a plane; 
here the critical point occurs at a = 0. We show this in Fig. Six curves forming three branches on Ba 
pass through a = at the angles given in eqs. ( 3.2.26| ). 

For g = 3, the locus Bu is shown in Fig. |l^. In this case, in addition to the eight phases that are 
contiguous at u = 0, there is a line segment on the real u axis from —0.573 to ~oo and a small O phase at 
the right end of this segment. There are also two complex conjugate pairs of O phases with the property 
that one pair is contiguous with the O2 and O2 phases while another pair is contiguous with the O3 , O3 and 
O4 phases. Normally, the fact that part of B extends to the point a = 0, as is the case here, means that 
for the given value of q the Potts antiferromagnet has a zero-temperature critical point. However, here one 
encounters noncommutativity in the definition of the free energy; if one sets a — (i.e., v = —1) first, then 
^L,5 vanishes and there is no such semi-infinite line segment on B. 

For (7 = 4, the locus Bu is shown in Fig. There is evidently another O phase that includes the negative 
real axis for u < —1. The line segment —1 < u < —0.5 on the negative real u axis and the small O phase 
at the right end of it are also present here. To show the behavior in the vicinity of the zero-temperature 
critical point of this g = 4 Potts antiferromagnet, we display the singular locus Ba in the a plane in Fig. |2^. 
The locus Ba passes vertically through the origin, a = and then curves into the Re(a) < half-plane. 




Figure 15: Locus Bui{L}) for the n — > 00 limit of the cyclic sqd strip with q = 10. Partition function zeros are shown for 
m = 20, so that Z is & polynomial of degree e = bm = 100 in v and hence, up to an overall factor, in u). 



3.5 Thermodynamics of the Potts Model on the Ly = 2 Strip of the sqd Lattice 
3.5.1 Ferromagnetic Case 



The Potts ferromagnet (with real q > 0) on an arbitrary graph has w > so, as is clear from eq. (1.5), the 
partition function satisfies the constraint of positivity. In contrast, the specific heat C is positive for the model 
on the (infinite-length limit of the) Ly = 2 sqd strip if and only if g > 1. For q ~ 1, fnq = 2.5 In a ~ 2.5K and 
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Figure 16: Locus l3u({L}): same as Fig. ^ for g = 2. 
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Figure 17: Locus Ba for the n — » oo limit of the 59,^ strip with q = 2. 
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Figure 18: Locus l3u({L}): same as Fig. ^ for g = 3. 
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Figure 19: Locus Bu({L}): same as Fig. ^ for 5 = 4. 
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Figure 20: Locus Ba for the n — > oo limit of the sq^ strip with g = 4. 



C vanishes identically. Since a negative specific heat is unphysical, we therefore restrict to real q > 1. For 



general q in this range, the reduced fi'ee energy is given for all temperatures by / = (1/2) In Ag^i as in ( 2.12 ). 
Recall that As,i = Xla- It is straightforward to obtain the internal energy U and specific heat from this free 
energy; since the expressions are somewhat complicated, we do not list them here. We show a plot of the 
specific heat (with fc^ = 1) in Fig. |^. One can observe that the value of the maximum is a monotonically 
increasing function of q. 

The high-temperature expansion of U is 

5J 
2^ 

For the specific heat we have 

5fcB(g- 1)A'2 



U 



1 



(9-1) 



K + OiK"^ 



2q2 



(5g + 14) 
5q 



K + 0{K^) 



The low-temperature expansions {K — > oo) are 



U = J 



+ 2{q-l)e 



-4K 



1 + 5e-^ + 7{q - 2)e-^^ + ©(e"^^ ) 



and 



,K^{q-l)e 



-4K 



1 
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-K 



+ !^(5_2)e-3^ + 0(e-4^) 



K 



as K 



(3.5.1) 



(3.5.2) 



(3.5.3) 



(3.5.4) 



4 4 

Comparing with our corresponding calculations for the {L^ — > oo limits of the) strips of the square and 
triangular lattices with the same Ly — 2 width, we can remark on some common features. In all of these 
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2.5 




Figure 21: Specific heat (with ks = 1) for the Potts ferromagnet on the infinite-length, width Ly = 2 strip of the sqd lattice, 
as a function of = J/{kBT). Going from bottom to top in order of the heights of the maxima, the curves are for q = 2,3, 4, 6. 
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cases, the high-temperature expansions have the leading forms 



[/ = -- 



AJ 



1 



(9-1; 



K + 0{K^) 



(3.5.5) 



where we recall that the coordination number is A = 3, 4 and 5 for these infinite-length strips of the square, 
triangular, and sqd lattices with width 2. (In the infinite- length limit, the longitudinal boundary conditions 
do not affect the coordination number.) Further, 



C 



2g2 



1 + OiK) 



For the low-temperature expansions for these strips. 



U = J 



-^ + 0((,-l)e-(^-i)^) 



and 



C (xkBK^{q-l)e 



-{A-1)K 



1 + 0(e 



as K —> 00 



as K ^ 00 . 



(3.5.6) 



(3.5.7) 



(3.5.8) 



In general, the ratio p of the largest subdominant to the dominant Aj's determines the asymptotic decay 
of the connected spin-spin correlation function and hence the correlation length 



1 

In p 



Since Xla and are the dominant and leading subdominant Aj's, respectively, we have 



Pfm 



X 



LA 



(3.5.9) 



(3.5.10) 



and hence for the ferromagnetic zero-temperature critical point we find that the correlation length diverges, 
as T ^ 0, as 

Cfm q^^e'^'^ as K ^ 00 (3.5.11) 

Comparing with the divergences in the correlation length at the ferromagnetic T ~ critical point that we 
have calculated for the infinite-length limits of the square and triangular strips with the same Ly — 2 width 
[p5|, nq, we see that all of these can be fit by the formula 



FAI 



-ie(^-i)^ as K 



00 



(3.5.12) 



3.5.2 Antiferromagnetic Case 

In this section we first restrict to the real range g > 4 and the additional integer values q = 2 (Ising 
case) and g = 3 where the Potts antiferromagnet exhibits physically acceptable behavior and then consider 
the remaining interval < g < 4, g 7^ 2, 3, where it exhibits unphysical properties. For g > 4, the free 
energy is given for all temperatures by ( 2.12| ), as in the ferromagnetic case but with J negative, and is 
the same independent of the different longitudinal boundary conditions, as is necessary for there to exist a 
thermodynamic limit. 

We show plots of the specific heat, for several values of q, for the Potts antiferromagnet on the (infinite- 
length limit of the) Ly — 2 strip of the sqd lattice in Fig In contrast to the ferromagnetic case, where the 
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positions of the maxima of C in the variable K increase monotonically as functions of q, the positions of the 
maxima of C in the antiferromagnetic do not have a monotonic dependence on q; they occur at K = —0.82 
for q — 2, K ^ —0.42 for q — 3 and K = —0.75 for q = 4, after which the values of the maxima occur at 
smaller values of with increasing integral q, e.g., K — —0.47 for q = 5 and K — —0.34 for q ^ 6). Recall 
that the antiferromagnetic Potts model with q = 2 and q = 3 involves frustration, and one can observe a 
rather different behavior in the specific heat for these values of q as contrasted with g > 4 in the plot. 




Fi gure 22: Specific lieat (with ks = 1) for tlie Potts antifcrromagnet on tlie infinite-lengtfi, widtli Ly = 2 strip of the 
sqd lattice, as a function of —K = —J/{kBT). Going downward in order of the heights of the maxima, the curves are for 
9 = 2,4,5,3,6. 



The high-temperature expansions of U and C are given by ( 3.5.l| ) and (3.5.2); more generally, these 
expansions also apply in the range < g < 4. As discussed above, the Ising case q — 2 is one of the 
cases where one must take account of noncommutativity in the definition of the free energy and hence of 
thermodynamic quantities. If one sets q = 2 first and then takes n oo, then / = /„g — (1/2) In Az,_4(g — 2) 
where \L,i{q = 2) was given in eq. ( ^.2.16| ), and the low-temperature expansions are 



U{q = 2) = - 



J 



1 + 126^^ + 36e''" + 80e^" + ©(e"^ ) 



5K 



and 



Ciq 



2„3K 



1 + 4e 



K 



100 



2K 



72 



as K ^ — oo 



as K 



(3.5.13) 



(3.5.14) 



9 9 

For g = 3, if one sets q — 3 first and then takes n oo, then / — fnq — (1/2) lnAL_4(g — 3), and the 
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low-temperature expansions are 



3^2 



K/2 



4e^ + ^e3^/2^0(e2^) 



as K 



and 



16 



3 4 



For the range <? > 4, the low-temperature expansions are given by 
(-J)e^ [,^^ 29g3 - 190^2 ^ 405(7 - 276 a- 



2(9-3)2 



(5(z-13) 



(9-3)2(q-2) 



and 



2„A' 



C 



kgK^e 
2(9-3)^ 



(5g-13) 



2(29g3 - 190g2 405q - 276) ^ 



(g-3)2(q-2) 



as K 



as iiT 



as X 



-oo 



-oo 



(3.5.15) 



(3.5.16) 



(3.5.17) 



(3.5.18) 



Note that for the antiferromagnetic case, U{T = 0) = for g > 4, but U{T = 0) = - J/2 = +\J\/2 for g = 2 
and U{T = 0) = -J/4 = -h| J|/4 for q = 3. The vanishing value of i7 at T = for g > 4 means that the 
Potts model can achieve its preferred ground state for this range of q, while the nonzero value of U{T = 0) 
for the Ising and q — 3 antiferromagnet is a consequence of the frustration that is present in this case. Note 
that the apparent divergences that occur as g ^ 2 and g 3 in eqs. ( |3.5.17 ) and ( 3.5.18 ) are not actually 
reached here since these expressions apply only in the region g > 4 (the discrete integral cases g = 2 and 
g = 3 were dealt with above). 

For the zero-temperature critical points in the g = 2 and g = 4 Potts antiferromagnet, 

Al,2 



PAFM,q=2,4 



X 



LA 



(3.5.19) 



Using the respective expansions ( p.2.24 )-( |3.2.25 ) and ( |3.2.27 )-(3.2.3C), we find that the correlation lengths 

(3.5.20) 



defined as in ( 3.5.9 ) diverges, as T 0, as 

£,AFM,q=2 



-3K 



and 



as K 



-oo 



lAFM,q=i ~ 1 ixti K ^ —OO . (3.5.21) 

Next, we consider the range of real < g < 4 aside from the integral case g = 2, 3. The first pathology 
is that the Potts antiferromagnet on the infinite-length limit of the cyclic Ly — 2 sqa strip has a phase 
transition at a finite temperature, call it Tp.^, while, in contrast, if one uses free boundary conditions, then 
either (i) there is no phase transition at any finite temperature, for 3<g<4orl<g<2, or (ii) there is a 
phase transition at a finite temperature Tp^s for 2<g<3or0<g<l, but Tpx 7^ 2p^5, so that there is no 
well-defined thermodynamic limit for the Potts antiferromagnet with non-integral g in the interval < g < 4. 
The Ising case q — 2 and g = 3 case have been dealt with in the preceding subsection. Concerning the value 
g = 1, as discussed earlier, one encounters nonconimutativity in defining the free energy. If one takes g = 1 
to start with and then n — > 00, the thermodynamic limit does exist, independent of boundary conditions, 
and / = fnq = 2.5K, U = — 2.5J = 2.5| J|, and C = 0. If one starts with g 7^ 1, takes n ^ 00, calculates 
fqn, and then takes g ^ 1, the thermodynamic limit does not exist since the result differs depending on 
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whether one uses free longitudinal boundary conditions or cyclic longitudinal boundary conditions. In the 
high-temperature phase, fqn = (l/2)lnAL,4, independent of longitudinal boundary conditions, but in the 
low-temperature phase, the expression for /g„ is different for the open and cyclic strips. There are also other 
unphysical properties, such as a negative specific heat and a negative partition function for certain ranges 
of temperature. 

Acknowledgment: The research of R. S. was supported in part at Stony Brook by the U. S. NSF grant 
PHY-97-22101 and at Brookhaven by the U.S. DOE contract DE-ACO2-98CH1O886.0 

4 Appendix 
4.1 General 

The Potts model partition function Z{G, q, v) is related to the Tutte polynomial T{G, x, y) as follows. The 
graph G has vertex set V and edge set E, denoted G = (V,E). A spanning subgraph G" is defined as a 
subgraph that has the same vertex set and a subset of the edge set: G" = (V^E') with E' C E. The Tutte 
polynomial of G, T{G,x,y), is then given by |^-[D 

r(G,x,y)= ^ (x-l)^(«')-^-(«)(y-l)^(«') (4.1.1) 

G'CG 

where k{G'), e(G'), and n{G') — n{G) denote the number of components, edges, and vertices of G', and 

c(G') = e(G') + fc(G') - n(G') (4.1.2) 

is the number of independent circuits in G' (sometimes called the co-rank of G'). Note that the first factor 
can also be written as {x — 1)''(<^)~''('5 ), where 

r{G) = n{G) - k{G) (4.1.3) 

is called the rank of G. The graphs G that we consider here are connected, so that k{G) — 1. Now let 

x^l + - (4.1.4) 

V 

and 

y = a = v + l (4.1.5) 

so that 

q={x~l){y-l) . (4.1.6) 

Then 

Z(G,q,i;)-(x-l)^-(«)(y-l)"(«)T(G,x,y) . (4.1.7) 
Note that the chromatic polynomial is a special case of the Tutte polynomial: 

P(G,(z) =(?^-(^)(-l)^-(^)+"(^)r(G,a; = l-g,y = 0) (4.1.8) 

(recall eq. (|l6|)). 

^Accordingly, the U.S. government retains a non-exclusive royalty-free license to publish or reproduce the published form of 
this contribution or to allow others to do so for U.S. government purposes. 
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Corresponding to the form ( 1.22| ) we find that the Tutte polynomial for recursively defined graphs com- 
prised of TO repetitions of some subgraph has the form 



Nx 



(4.1.9) 



4.2 sQd Strip with Free Longitudinal Boundary Conditions 

The generating function representation for the Tutte polynomial for the open strip of the sqd lattice comprised 
of TO + 1 squares with edges joining the lower-left to upper-right vertices and the upper-left to lower-right 
vertices of each square, denoted Sm, is 



We have 



where 



with 



and 



with 



where 



TT{Sra,x,y;z) = ^ T{Syn,x,y)z'^ . 

Tt[S, X, v; z) = — —-^ 

VT{S,x,y]z) 

Mt{S, X, y; z) = At,s,o + At,s,iz 

Ar^sfi = 2x + 4a;?; + 3x^ + 2y + 3y^ + y^ + x^ 
^T,s,i — xy{y + y'^ + X — 2xy^ + xy + x^ — 2yx^ — x^y'^) 



T^T{S,x,y;z) = - ^T.s,jz) 



•^T,S,(l,2) — 2 



r.ix + 3) + + 2?/ + 3) + 2 ± /R^ 



(4.2.1) 

(4.2.2) 

(4.2.3) 

(4.2.4) 
(4.2.5) 

(4.2.6) 
(4.2.7) 



Rt = 4 + 12a; + 12y + 22xy + lix^ + 21y^ + &x^ + 20y^ + IGxy^ 



+\Ox^y + X* + IQy" ~ Ax^y^ - 2y'x + Ay" - 2x^y^ + / 



(4.2.8) 



The corresponding closed- form expression is given by the general formula from |22| , as applied to Tutte, 
rather than chromatic, polynomials, namely 



T{Sm,x,y) = 



S,2 



(At,s,i)™ + 



T,S,2 — At,S,1 



(At,S,2)'^ 



(4.2.9) 



It is easily checked that this is a symmetric function of the Xsj, j — 1, 2. 
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4.3 Cyclic and Mobius Strips 

We write the Tutte polynomials for the cyclic and Mobius strips of the sqij lattice with width Ly = 2 as 

5 

T{L^,x,y) = ^ctxj(At,l,,)" (4.3.1) 

where it is convenient to extract a common factor from the coefficients: 

CT,.. ^ . (4.3.2) 

Of course, although the individual terms contributing to the Tutte polynomial arc thus rational functions of 
X rather than polynomials in x, the full Tutte polynomial is a polynomial in both x and y. We have 

At,l,i = 2 (4.3.3) 

and, with 

7t,23 = 2/' + 2i/ + 3y + 2x + i (4.3.4) 
i?r,23 = 16 + 16a; + 32y + 12xy + 33y^ - ixy^ + lOy* + ix^ + 20y^ + iy^ + y^ (4.3.5) 

the results 



A- 



T,L,(2,3) 



23 



(4.3.6) 



and 



The coefficients are 



Ar,L,(4,5) = At,s,(i,2) (4.3.7) 

CT,L,i = i(x - l)(y - \){xy -x~y-2) (4.3.8) 
ct,l,2 = ct,l,3 ^xy-x~y (4.3.9) 

CT,LA = CT,L,5 = 1 . (4.3.10) 

These are symmetric under interchange of a; <-> y, which is a consequence of the fact that the c^j are 
functions only of q. 

4.4 Special Values of Tutte Polynomials for Strips of the sqd Lattice 

For a given graph G = (V,E), at certain special values of the arguments x and y, the Tutte polynomial 
T{G,x,y) yields quantities of basic graph-theoretic interest §-||l^, |3^. We recall some definitions: a 
spanning subgraph G" — {V, E') of G is a graph with the same vertex set V and a subset of the edge set, 
E' C E. Furthermore, a tree is a graph with no cycles, and a forest is a graph containing one or more trees. 
Then the number of spanning trees of G, Nst{G), is 

7Vst(G) =r(G, 1,1) , (4.4.1) 

the number of spanning forests of G, Nsf{G), is 

7Vsf(G) =r(G, 2,1) , (4.4.2) 
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the number of connected spanning subgraphs of G, Ncssg{G), is 



Ncssg{G)=T{G,1,2) , 
and the number of spanning subgraphs of G, Nssg{G), is 

^ssg(G) =T(G,2,2) . 
From our calculations of Tuttc polynomials, we find that 

NsASm) = 22("+2) • 3" 



(4.4.3) 



(4.4.4) 



(4.4.5) 



NsFiSm) = (19 + liVS )(9 + 5\/3 )™ + (19 - llV3 )(9 - bVs f' 



(4.4.6) 



NcssaiSrn) = ( 19 + ) [2(7 + VA6 )]" + (l9 - ) [2(7 - V46 )Y 



23 



23 



;)5m+6 



(4.4.7) 



(4.4.8) 



For the cyclic Ly = 2 strip of the sq^ lattice, Lm, we first note that for m > 3, Lm is a (proper) graph, 
but for m = 1 and m = 2, L,„ is not a proper graph, but instead, is a multigraph, with multiple edges (and, 
for m — 1, loops). The following formulas apply for all m > 1: 



iVsT(im) = 22("-l) . 3" 



TO 



(4.4.9) 



(9 + 5V3 y" + (9 - 5V3 y" - (^7 + 3\/5 (7-3V5 



(4.4.10) 



NcSSciLrn) = S ■ 2™"! 

7V46 



-| m r 

2(7 + \/46) + 2(7-Vl6 



TO 
2" 



m — 1 



7V46 
23 



2( 7- \/46 



m — 1 



^5m 



This result, eq. ( 4.4.12 ), is a special case of a more general elementary theorem, namely 

Nssg{G) - 2^(G) 



(4.4.11) 



(4.4.12) 



(4.4.13) 



This is proved by noting that a spanning subgraph G' G consists of the same vertex set F as G and a 
subset of the edge set E of G. One can enumerate all such spanning subgraphs as follows: for each edge in 
E, one has the option of including or excluding it, keeping the other edges fixed. This is a twofold choice 



for each edge, and the result (4.4.13) therefore holds. This general result subsumes the previous specific 
relations NssG — 2^™ and NssG = 2^™ for the cyclic strips of length — m and width Lj, = 2 of the 
square and triangular lattices. We recall that the number of edges is given by e(G) = (i) 3to, (ii) 4to, and 
(iii) 5to for the — Ly — 2 strips of the (i) square (ii) triangular, (iii) sqd strips, respectively, while the 
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number of vertices is given by n = v{G) = 2m for all of these strips. Let us denote re = e{G)/n{G), whence 
re{G) = (i) 3/2, (ii) 2, (iii) 5/2 and NssciLc^m) = 2'-'=('5'" for these strips. 

Since T{G„nX,y) grows exponentially as m — > oo for the families Gm = Sm and Lm for {x,y) ~ (1, 1), 
(2,1), (1,2), and (2,2), one defines the corresponding constants 

Zset{{G})^ lim n{GyHnN,et{G) , set ^ ST, SF, CSSG, SSG (4.4.14) 

n{G) — ^oo 

where, as above, the symbol {G} denotes the limit of the graph family G as n{G) — > oo (and the z here 
should not be confused with the auxiliary expansion variable in the generating function (4.2.1) or the Potts 
partition function Z{G,q,v).) General inequalities for these were given in |l5[] . 
Our results yield 

^st({G}) = ln2 + iln3 ~ 1.242453.. for G = S,L (4.4.15) 



ZSF 



({G}) = iln(9 + 5V3 ) ~ 1.435658 for G = S,L (4.4.16) 



2 



zcssg{{G}) ^ i ln[2(7 + \/46 )] ~ 1.658267 for G = S,L (4.4.17) 

and 

zssg{{G}) = ^ In 2 ~ 1.732868 for G = S,L (4.4.18) 

In Table I we summarize the results for Zs{{G}) for the infinite-length limit of the width Ly ~ 2 strips of 
the square, triangular, and sqd strips (which are independent of the longitudinal boundary conditions). In 
this table we include a comparison of the exact values of zst that we have calculated with the upper bound 
(u.b.) for a /c-regular graph |Q 

ZST < zsT,u.,. , ZST..U.,. = In ( [fc(f_'2)p/2)-i ) (4-4-19) 

The cyclic Ly — 2 strips of the square, triangular, and sqd lattices are fc-regular, with /c = 3, 4, 5, respectively. 
One thus has 

zsT^u.b. {sq{Ly = 2))= ln(4/V3) (4.4.20) 

zsT.u.b. {tri{Ly = 2)) = 3 ln(3/2) (4.4.21) 

zsT,u.b.{sqd{Ly = 2)) = 81n2- ^lnl5 (4.4.22) 

For this table we define 

rsT{{G}) = . (4.4.23) 

As is evident from the table, zst{{G}), zsf{{G}), and zcssg{{G}) increase as one goes from the Ly — 2 
strip of the square strip to that of the triangular, and sqd strip, as the degree increases from 3 to 4 to 5. 
This also follows from (4.4.13) for zssg- A similar dependence on degree (coordination number) was found 
for ZST in 11,0. 
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Table 1 : Zs ({G}) for the — ► oo limit of the width Ly = 2 strips of the (i) square lattice, (ii) triangular lattice (constructed 
by starting with a square strip and adding edges joining the lower left to upper right vertex of each square, and (iii) sg<j lattice, 
i.e., square lattice with next-nearest- neighbor couplings. 



^s{{G}) 


{G} = sq 


{G} = tri 


{G} = sqd 


zst{{G}) 


(l/2)ln(2 + V3 ) 
= 0.658479 


(l/2)ln[(7 + 3V5)/2] 
= 0.962424 


lii(2\/3 ) 

= 1.242453 


rsT{{G}) 


0.7867 


0.7912 


0.8377 


zsf{{G}) 


(l/2)ln[2(2 + V3 )] 
1.005053 


(l/2)ln[2(3 + 2V2 )] 
= 1.227947 


(l/2)ln(9-F5\/3 ) 
= 1.435658 


zcssg{{G}) 


(l/2)ln[(5 + Vl7)/2] 
= 0.758832 


(l/2)ln[2(3 + 2V2 )] 
= 1.227947 


(l/2)ln[2(7 + v'46 )] 
= 1.658267 


zssg{{G}) 


(3/2) ln2 
1.039721 


21n2 
= 1.386294 


(5/2) In 2 
= 1.732868 
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